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PHOTON COUNTING: A PROBLEM IN 

CLASSICAL NOISE THEORY 

By Sherman Karp and John R. Clark 
Electronics Research Center 

SUMMARY 

In this report we formulate the general problem of determining 
the photoelectron "counting" distribution resulting from an elec- 
tromagnetic field impinging on a quantum detector. Although the 
detector model used was derived quantum mechanically, our treat- 
ment is wholly classical and includes all results known to-date. 
This combination is commonly referred to as the semiclassical 
approach. The emphasis, however, lies in directing the problem 
towards optical communication. 

The electromagnetic field is assumed to be the sum of a 
deterministic signal and a zero-mean, narrow-band, gaussian random 
process, and is expanded in a Karhunen-Lodve series of orthogonal 
functions. Several examples are given. It is shown that all the 
results obtainable can be written explicitly in terms of the noise 
covariance function. Particular attention is given to the case 
of a signal plus white gaussian noise, both of which are band- 
limited to ± B Hz. Since the result is a fundamental one, to 
add some physical insight, we show four methods by which it can 
be obtained. Various limiting forms of this distribution are 
derived, including the necessary conditions for those commonly 
accepted. The likelihood functional is established and is shown 
to be the product of Laguerre polynomials . For the problem of 
continuous estimation, the Fisher information kernel is derived 
and an important limiting form is obtained. The MAP (maximum a 
posteriori) and ML (maximum likelihood) estimation equations are 
also derived. In the latter case the results are also functions 
of Laguerre polynomials. 


INTRODUCTION 


Since the advent of the laser, a problem of growing impor- 
tance in optical communications and coherence theory has been 
the determiniation of the output statistics of a quantum detector 
excited by a narrow-band gaussian source. In optical communica- 
tions, knowledge of these statistics is necessary for the appli- 
cation of the techniques of optimum detection and estimation 
theory. In the physical theory of coherence, these statistics 
are a means by which the light incident on the detector can be 
studied. In both cases a useful statistic, which is relatively 
easy to evaluate, is the probability p N ^_(k) of detecting k events, 
or "counts," in the time interval (0,t]. In an idealized detector 
the conditional probability of k counts in (0,t], given the in- 
cident radiation, can be shown to obey a Poisson law, with the 
time-averaged intensity of the field as rate parameter. 

Although this idealized detector model is based on the quantum 
theory of photodetection, statistics such as p N ^_ (k) can, if de- 
sired, be correctly calculated by using only classical tools, 
provided the proper quantum model has been postulated. This is 
the essence of the so-called "semiclassical" approach, which we 
shall use in this report (ref. 1) . The value of this approach 
lies not so much in its freedom from quantum-mechanical subtleties, 
as much as in the ease with which it allows meaningful physical 
interpretations to be made and comparisons to be drawn. We shall 
therefore not dwell on the physical considerations leading to the 
derivation of the Poisson model, since it is generally accepted 
as an adequate probabilistic description of an ideal quantum 
detector localized at a point in space (ref. 2) . The physics 
underlying this model and others is discussed at length in the 
literature on quantum detectors. 

The problem of determining the counting distribution p N ^(k) 
has been approached by a number of authors, within the framework 
of laser intensity fluctuation studies (refs. 3 through 8). It 
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is intimately connected with the classic noise theory problem of 
finding the statistics of the time-averaged square of a real 
process. Random functionals of this type have been treated in 
considerable generality for gaussian processes; we list only a 
few of the more important works in references 9 through 15. 

Here, we present a general method for finding p N ^(k) when 
the light incident on the detector is a deterministic signal 
plus a narrow-band gaussian process. Our results are shown to 
encompass as special cases several previously known results. 

The method is a generalization of Mandel's (ref. 5) in which we 
allow the incident field to have a nonzero mean, and p N ^_(k) is 
expressed in terms of the cumulants of the time-averaged intensity. 
Several specific examples are worked out in detail, which are 
in agreement with known results when such results exist. Con- 
siderable attention is given to the important case in which the 
radiation is a deterministic signal in band-limited white gaussian 
noise. For this example, we show how the results can be applied 
to some problems of interest in optical communications . 

DERIVATION OF THE COUNTING DISTRIBUTION 

First, we consider the Poisson model for the ideal quantum 
detector. Given a counting statistic N^_ at time t, conditioned 
on a function [a(x); 0 < x < t] , such that it obeys a Poisson 
law. 


C 


Pr{N t = k| [a(x);0 < x 5 t] } 


_ [m(t)] k „-m(t) 
“ ' k! 


(1) 


m (t) 


a 


f~ | a ( x ) 
0 


2 . 
dx 


(2) 


we want to know the probability Pr{N fc = k} when a(x) is the complex 
envelope of a narrow-band gaussian process. Here a = n/hv, where 
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n is the quantum efficiency, h is Planck's constant, and v is the 
frequency. With m^ the random value assumed by m(t') at time 
t' = t, it is clear from Eq. (1) that the probability of k counts 
in (0,t] is given formally by 


Pn t (k) = in- E { e t ”t} (3) 

where the expectation is taken over the random variable m t » 

We shall take a(x) to be the complex envelope of a real, not 
necessarily stationary, gaussian process a(x), 

a(x) = Re[a (x)e^ 27rf o T ] (4) 

which is assumed to be narrow-band about some very high frequency 

f . We shall also assume that the covariance function of a(x) is 
o 

real. With these assumptions we can show that the real and 
imaginary parts of a(x) are independent and have the same covariance 
function. In addition we assume that a(x) can be written 


a (x) = s (x) + n (x) 


(5) 


where s(x) is a deterministic signal and n(x) is a zero-mean 
gaussian random process. 

Instead of evaluating Eq. (3) directly for p Nt (k), we shall 
find the characteristic function of N t : 


M (jv) 
t 



( 6 ) 


{ m (e jv - 1) 
Eje 


(7) 


Then p Nt (k) can be found either by direct inversion of M Nt (jv), 
or from one of the following formulas: 


1 


4 





5 = 1 


(8) 



(k) 


k! 




Un(l - ?)] 




(-l) n 

n! 


^k+n 


(9) 


vu is the i-th moment of m t ; in terms of M Nt (jv), 

y = 1^- M [An(l + 5) ] (10) 

1 9? 1 w t 5=0 

Equation (8) is the equivalent of differentiating the probability 
generating function of N^., and Eq. (9) is the result of expanding 
the exponential in Eq. (3) in a power series. 

We see from Eq. (7) that f%^_(jv) is simply the moment 
generating function of m^, 


M (u) 
m t 



( 11 ) 


evaluated at u = e-* v - 1; thus, we can confine our attention to 
the random variable, m^. 

It is convenient to expand ad) in a Karhunen-Lo^ve series 
on [0 , t] (ref . 16 ) : 






a (t) = 2 a i d) 
i=l 1 1 


= 53 (s i + n ± ) 4> i (t) 


( 12 ) 


iL 


The equality, of course, is in the sense of "limit-in-the-mean. " 
The coefficients are given by 
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(13) 


n ± = (n, <J) i ) 


s i = (s, <j> i ) 


(14) 


where 


(x, y) = f x (t) y* (t)cIt 

* / 0 


(15) 


and the { > are eigenfunctions of the integral equation 


'i+i< u > - Vi -/ 


K„(u, 


v) <J> i (v) dv 


(16) 


Here K (u, v) = E{n(u)n*(v)} is the covariance kernel of the noise. 
The {<|k} are normalized so that (<Jk, 4> j ) = • 

It is clear from the orthonormality of the eigenfunctions 
that m t can be written 

m , = a 23 | | 2 = ot 23 I s i + n i I ^ (1^) 

z i i 


Notice that m t is the energy in the process at time t. Since n(x) 
is a zero-mean gaussian random process, the {n^ are gaussian 
random variables, with E(n^) = 0 and E (n^n^ ) = A ^ 6 ^ j . This 
orthogonality depends critically upon choosing the basis to 
satisfy Eq. (16) uniquely. If, however, there is no noise 
(n (x ) = 0), a(t) = s(t) can be expanded in any complete ortho- 
normal set {iJk} on [0 , t] , and Eq . (16) is irrelevant. For this 

case, with c ± = (a, \p ± ) , Eq. (1) can be written f 


Pr 


|N t = k | [a (x) ; 0 < x < t] J 

[e a| Ci | 2 j k 


exp - 23 a I c i 


(18) 
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p 


Each coordinate axis in the space contributes an independent 
Poisson variate , 


Pr { N ti = k j_U a ( T )'* 0 < T 1 t]j 



(19) 


where N. = S N . and k = 2 k. (ref. 17). This is clearly 

L. 2 _ tl 1 1 

independent of the particular basis chosen; only the {c^} change. 

Each axis always contributes an independent Poisson variate. In 

addition, since S |c. is the energy of a(x) in [0,t] , the 

i 

conditional density of is independent of the functional form 

of a(x). If we choose {^.} to be the sinusoidal set on [0,t], 

2 

then | c^ | represents the energy of a(x) at the frequency of 
If a(x) = ipj (x) for some j, then k = k^ and Eq. (18) reduces to 
Eq. (19) with i = j. 

For the more general case of nonzero noise, we observe that 
only one particular orthonormal set can be used as a basis for 
expanding a(x), if we desire each axis in the space to contribute 
an independent variate to N.. That basis, of course, must 
satisfy K n 4> = X<f), Eq. (16). (The only exception is white gaussian 
noise, for which Eq. (16) is satisfied by any complete orthonormal 
set; then each axis in the space contributes an independent, 
identically distributed random variable.) 

For narrow-band gaussian noise, we now show that the con- 
tribution from each axis is an independent Laguerre-distributed 
variate and, consequently, that can always be represented as 
the sum of independent Laguerre random variables. 

The real and imaginary parts of n.^ = (n, cf>.) are independent, 

each with variance A^/2, thus (ref. 18, p. 196), 
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2 2 
S ± | + A 


P|a.| (A) = < 


2A 

X i 


1 T °( 2 ^r A )' 


A k 0 


0, elsewhere 


( 20 ) 


v.. 


the so-called Rician density, where I is the zero-order modified 

® 2 

Bessel function of the first kind. It follows that | a^ | has the 
density function 

I s ± | 2 + A 


P, |2 (A) =< 


k : ‘‘ 


A > 0 


0, elsewhere 


( 21 ) 


Since the {a^} are independent, it follows that the {|a^|^} are 
also independent. Thus, from Eqs. (11) and (17) we see that 


M (u) = II E 
m t i 


a|a i | u 


( 22 ) 


Now, with x a real, gaussian variate, E|e X u | is given by (ref. 
19, p. 396) 


E(e 


2 

x u 


exp 

E 2 (x) U 

.1 - 2var(x)u_ 

1 1 - 2var(x)u 

1/2 


^Re u < 2var(x)^ 


(23) 


Thus , 


M m t < u > 


= n — a A : ~ u ex P 
1 


a| Si | 2 u / 1 \ 

t t , ( Re u < t — I (24) 

1 - aA . u ' 1 a max X . I 

1 V i 1 / 
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and 


as 


m n (jv) = n 
t i 


(e^ v - 1) 


1 - aX^ (e- ,v - 1) 


exp 


1 - aX ± (e :5V - 1) 


(25) 


The quantity within the brackets is the characteristic 
function of a Laguerre random variable for which the probability 
distribution is (ref. 20) 



(26) 


Lj, is the zero-order Laguerre polynomial (ref. 21, Eq. (8.97)). 
Thus, we see that N t can be represented as the sum ^ of 

Laguerre variates, as asserted. If the noise is gaussian, thej-e 
always exists an operation which acts on Eq. (18) , producing 
another independent coordiniate system of orthogonal axes 
(Eqs . (25) and (26)). 

If s. = 0 for some i. 



Hence any coordinate <f>^ for which s^ = 0 contributes a Bose- 

Einstein random variable to N . . 

t 

Proceeding formally, we note that p N ^_(k) is just the in- 
finite convolution of the {p M (k • ) } for all i. As this is 

LN ti 1 

usually difficult to evaluate explicitly by using Eq. (26) , we 
shall concentrate our attention on Eq. (24) instead. 
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We mention in passing that if the noise n(x) is not gaussian, 
the interpretation in the preceding three paragraphs cannot be 
used; i.e., there does not exist an orthonormal basis in which 
a(x) can be expanded such that N t can be represented as the sum 
of independent variates. 

M mt (u), Eq. (24), can be cast into a completely equivalent 
form that does not explicitly involve the eigenvalues { A^ } (refs. 

9 and 22) . The identities which make this possible are (ref. 23) 



TrK 


(k) 


n 




(f , 


t’> df 


and 


(28) 


£ 


2 A k 

l 


(s, 


K (k) s) 

n 


(29) 


where 


K^ k) (u, v) = J K^ k 1) (u,?)K n (?, v) d£ 


(30) 


K (1) = K 
n n 


= identity operator 


(31) 


(32) 


-1 . 


By expanding log (1 - aA^u) and (1 - aA^u) in power series, and 
using Eqs . (28) and (29), we get the immediate result 


M (u) = exp £ 
m t k= 


. K k k / l \ 

-k T U ' U " a max (A T)J 


(33) 


and {k^}, the cumulants of m^, are 

K ± = j\i - 1) :TrK j J l) + i! (s,K^ 1_1) s)J a 1 


( 34 ) 
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J 


For the case, s(x) = 0, a number of authors have obtained 
expressions equivalent to Eq. (24) or Eq. (33) (refs. 3, 4, 7, 
and 9) . Related results have also been obtained for real, rather 
than complex, a(x) (refs. 9, 10, 11, 12, 14, and 15). 

The form of the noise will often dictate which of the two 
formulas, Eq. (24) or Eq. (33), will be the more convenient in 
practice. The integral equation (Eq. (16)) is difficult to 
solve, and even when the eigenvalues are available, Eq. (24) 
might not simplify significantly. On the other hand, as Slepian 
(ref. 11) has pointed out, the iterated kernels (Eq . (30)) are 

usually difficult to evaluate beyond the first few orders; 
however, they are often all that is needed for a good approxima- 
tion. 

We can now find p^^(k) by using one of the formulas, Eq. 

(8) or Eq. (9). It turns out that Eq. (9) can be evaluated 
almost trivially with the help of the moment generating function 
M mt (u) . Since the moments {y^} and the cumulants {k^} are 
related, we need not evaluate Eq. (10) directly. We use instead 
the formula (ref. 24), 



k. > 1 

3 


Along with Eqs . (34) and (35), Eq. (9) gives an exact expression 

for P£j t (k), the desired counting distribution, in terms of the 
covariance K n and the signal s. 

I The fundamental quantities here are the cumulants {k^}, in 

terms of which a number of the statistical characteristics of 
N^_ can be expressed. Some of the more important relations are 
derived in the Appendix. Using Eqs. (28), (29), (34), and (A2) , 

we find the mean and variance of to be 


11 


(36) 


= aTrK n + a(s, s) 


var (N t ) 


23 + a(s, s) + £ [ a *jJ 

+ 2 23 [ aX i] ® I ® I 


( 2 ) 


= aTrK fl + a ( s , s ) + a TrK n 


2a 2 (s , 


K s) 
n 


(37) 


Notice that the only quantities involved in Eqs. (36) and (37) 
are the noise covariance, the signal, and the time, t. 

It should be emphasized that our method of finding p^ t (k) 
through the characteristic function is not always the easiest 
path to the desired answer. Indeed, for certain noise sources 
the most direct route to PN^fk) might be the direct evaluation 
of P mt (M) , the density function of m t , and use of Eq. (3). 
Whichever method is to be used is best decided when the precise 
form of K n has been ascertained. 

EXAMPLES AND APPLICATIONS 


It can be seen that the mathematical form of p N ^_(k) is 
considerably more complicated than that of (jv) — not a 
surprising situation in nonlinear noise problems. In most 
applications of interest, one expects that p N (k) will be much 
easier to find by evaluating M N (jv) first, rather than by going 


directly to Eqs. (9), (34), and (35) 

follow, this is certainly the case. 


In the examples to 



Noise Only (s(t) = 0) 

2 

Integrated white noise: K n (u, v) = p min (u, v).- Although 

this example appears to be only of academic value, it is illus- 
trative of a particularly simple means of obtaining p Nt (k). It 
is well known that the eigenvalues of Eq. (16) are (ref. 19, 
page 196) 


X i = 


4p 2 t 2 


TT 2 (2i-l) 2 ' 


i = 1, 2, 3, . . . 


(38) 


Inserting this in Eq. (24) and using a tabulated product rule, 
we get (ref. 21, Eq. (1.4313)) 


M (u) = sec (pt/au) (39) 

t 

The coefficients in the Taylor series expansion of Eq. (39) are 
the moments {y^} of m^.; since (ref. 21, Eq. (1.4119)) 


sec x 


E 


k=0 


l E 2kl ,2k 

(2k): 



we have 


(40) 


y 


k 



(ap 2 t 2 ) 


(41) 


and the {E.^} are the Euler numbers. Equation (9) now gives p Nt (k) 
directly: 


N , 


(k) 


oo 


E 

n=0 


(~l) n 
(2k+2n) ! 



E 2k+2n 


f 2 2 x k+n 
(ap t ) 


(42) 


Because of the simple form of M^fu), we can find the 
counting distribution without first having to compute the curnu- 
lants {k ± }. The mean and variance of the counts work out to be 


13 


(43) 


i a 2,2 
E (N t ) = 2 P t 


, a 2,2 , a 4 4 
var(N fc ) = ■j P t + g— p t 


(44) 


First-order Markov noise: K n (u,v) 


= P exp - 3 I u-v | . - 


Equation (24) becomes 


M (u) 
m t 



1 

aA^u ' 


(45) 


with A^ a solution of a transcendental equation (ref. 9). 
Equation (45) has been evaluated by an indirect technique (ref. 
14) and is given by 


M m (u) 
m t 



+ 1 1 Ip/ - sinh [4 ' TT u )]} 1 <46) 

This does not easily yield a useful expression for p N ^_ (k) ; 
however, using the Leibnitz differentiation rule, Bedard (ref. 3) 
has obtained recurrence relations for the counting distribution 
and its factorial moments. For certain limiting cases, approxi- 
mate counting distributions have been obtained (ref. 7). 

The mean and variance of N t are easily found to be 

E(N fc ) = aPt (47) 


var(N^_) = aPt + ^23t + e - lj (48) 

Band-limited white gaussian noise .- Let n(x) be a white 
gaussian process, band-limited to ±B Hz, with two-sided spectral 
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height N q . Taking the first 2Bt + 1 eigenvalues of Eq. (16) to 
be the same (N q ) , and the rest to be zero (ref. 19, page 193), 
Eq. (24) becomes 


M 


that is. 


v u) = 


2Bt + 1 


(49) 


V jV) = [l - aK . J 


2Bt + 1 


(50) 


This can be recognized as the characteristic function of the 
negative binomial distribution (ref. 17) 


/2Bt + k\/ , \ 2Bt + 1 / aN \ K 

PN t°° \ k /V + aN o) V 1 + aN o) 

If 2Bt << 1 (t << -^g) , only one eigenvalue of Eq, 
significant, and Eq. (50) reduces to 


(51) 

(16) is 


M (jv) 
t 


1 

1 - aN Q (e^ V - 1) 


(52) 


In other words, the counting distribution is Bose-Einstein : 

/ aN \ k 

P N t ^ = 1 + aN Q (l + aN Q j (53 ^ 

in agreement with Eq. (27) . For 2Bt >> 1, and aN << 1, it can 

° . * 

easily be shown that p N ^(k) approaches a Poisson distribution. 

On the other hand, for a nonzero signal and no noise, p N ^_(k) is 
also Poisson, as we have seen. This behavior, in the case of 
band-limited white-gaussian noise, is due to the complete 

Note: This limiting form can also be shown to be true for first- 
order Markov noise, Eq. (46) , when Bt >> 1 and B >> 4Pa. 
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"smoothing out" of the intensity fluctuations, while in the 
signal only case, it is due to the deterministic nature of the 
signal. 

Signal Plus Noise 

Now we remove the restriction that s(t) be zero. Assume 
that n(x) is band-limited white gaussian noise, as above. In 
addition, assume that s(t) is band-limited to ±b' Hz, b' < B. 
Equation (25) applies with - N q/ and P^ t ( k ) is just the 
(2Bt + l)-fold convolution of the Laguerre density (Eq. (26)) 
with itself. Omitting the details, we note that the result is 
(ref. 21, Eq. (8.9771)) 


N, 


(k) 


(aN Q ) 


(1 + aN ) 
o 


k + 2Bt + 1 


a (s , s ) 
1 + aN 

o 




(s, s) I 

(1 + aN Q )J 


(54) 

2Bt 

where L^. is the 2Bt-order Laguerre polynomial. The mean and 
variance of N^. are found to be 


E(N t ) = (2Bt + l)aN Q + a (s , s) 


(55) 


var(N t ) = (2Bt + 1) (1 + aN Q )aN o + (1 + 2aN Q )a(s, s) (56) 


Equation (54) could have been found by evaluating p m (M) 


first, and using Eq. (3) — the alternative approach suggested 
earlier. For this example, (M) turns out to be a "noncentral 
chi-square" density (i.e., the density of the sum of 2Bt + 1 in- 
dependent Ricean variates) : 


p m (M) 
m t 


I 1 [Ml 

= <N a L(s, s ) aj 


M + (s , s ) a 
Bt aN 

e ° I 


( 2 


2(s, s) 1/2 1/2 

2Bt l 1727; M 
a ' N^ 
o 


) 


0, elsewhere 


M > 0 


(57) 
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p Nt (k), Eq. (54), results when the expression 



1 

kJ 


M k e" M 


p (M) dM 
m t 


( 58 ) 


is evaluated. 


For 2Bt << 1, Eq. (54) reduces to Eq. (26); for 2Bt >> 1, 
however, PN t (k) approaches a Poisson distribution when aN Q is 
small. This is clear from the following discussion. 

Let aN Q << 1 and 2Bt >> 1; then, asymptotically. 



-2BtaN 

o 

e 


(59) 


a(s, s) 
1 + aN 

o 


a (s , s) 


aN 

o 

1 + aN 

o 


aN 

o 


2Bt 


[- 


(s, s) 


V 1 


+ aN 


o>] ~ ( : 


2Bt + 


(s , s) 
N_ 




(60) 


(61) 


(62) 


These Equations, together with Eq. (54), yield 

fa(s, s) + 2BtaN l k -fa(s, s) + 2BtaN 3 
~ k! ° (63) 

Thus for small aN Q and large 2Bt, the counting distribution is 
Poisson, with the rate parameter the sum of a signal intensity 
and an independent noise intensity. Expressing this in another 

way, we note that N. is the sum of two independent Poisson 

*- 1 
variates; one associated with the signal (intensity ^(s, s)a). 
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and one associated with the noise (intensity 2BN a) . 

o 

Our results for band-limited gaussian noise are well known 
in the field of coherence theory (refs. 1, 2, 5, 6, 8, 20, and 
25); in optical communications, however, they are just beginning 
to find application (refs. 26, 27, and 28). In earlier papers 
the asymptotic form (Eq. (63)) was used, with heuristic justi- 
fication (refs. 29 and 30). 

ESTIMATION FOR CONTINUOUS WAVEFORMS 

For band-limited white gaussian noise, simple detection 
and estimation problems can be solved with relative ease because 
of the independence of noise samples taken at the Nyquist rate. 
Assume that the interval [0,t] is broken into 2Bt equal sub- 
intervals, each of length 1/2B sec. Then the probability of k^ 
detector counts in the i-th subinterval is given by 


k. a ' S il 

(ay 1 1 + aN o r Is. I 2 1 

1 + aN J 1 + k i k i L N o (1 + a yJ 


(64) 


where | s^J 2 = | s I ^/2B. Moreover, the counts in different 
subintervals are independent, because the noise samples are 
independent. Thus the joint probability p^ (k) of the counts in 
each subinterval is simply the product of probabilities (Eq. (64)) 


P N < k > = 


2Bt 

n 

i=l 


k. a|s il 

(otN o ) 1 " r+~5S£ ^ [ Is. | 2 ] 

II + aN ) X + k i e Lk i L V 1 + aN o } J 


( 65 ) 


with N = vector of 2Bt counting observables; 
k (k^, ^2 7 • • • r k 2Bt ^ * 
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It is important to note that, for this representation to be 
valid when the signal waveform is band-limited to ±B Hz, it is 
essential that the subintervals be of length precisely 1/2B sec. 

Equation (65) suggests that we define the likelihood func- 
tion A (k f s_) to be 


A(k|s) 


2Bt 
n e 
i=l 


a|s i 


2 


1 + aN 

o 



l3il 1 

■ N o ~(l + aN 0 )_ 


( 66 ) 


It is now straightforward to set up decision structures for the 
multiple-hypothesis detection problem. This has been done 
elsewhere for a number of signaling schemes (ref. 27), so we 
shall restrict our attention to the estimation problem. 

We start with the assumption that s is a transformation, 
s[t, x(t)], of some function, x(x), which we wish to estimate; 
x can be a modulation, a set of parameters, or a single number. 
As is well known, a necessary condition for x^^ (t) to be the 
maximum a posteriori (MAP) estimator of x(x), 0 < x < t, is 


V 

x 


f 


SLnA + &np 


x (x) 


[X (x ) ] 

J x 


X MAP 


0 


(67) 


where V x is the gradient operator with respect to x(x), and 
P x ( T ) t x (x) ] is the a priori probability density functional of 
x(x). In the absence of prior information, of course, Eq. (67) 
reduces to the maximum likelihood equation. 


V x £nA 



( 68 ) 


With some manipulation, V InA works out to be 

X 
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V SLnA 
x 


2a 


1 + 


2Bt 

SHT £, c lk 

O 1=1 


. , s . ) Re 

l ' l 


[> _!=il 

[fi 3 x(t)J 


(0 < T £ t) 


where 


c(k i f s i ) 




N o (l + 1 aH 0 )j 

ISil 2 "-. 
V 1 + aN 0 )J 


1 


( 69 ) 


(70) 


The problem is now simplified if we consider x as a set of para- 
meters (say, M of them). Then, with x = (x^, , . .., x^) , 

Eq. (69) reduces to the set of M expressions: 


V ink 
x 


2a 


1 + aN 


2Bt 

2 c(k, , 


o i=l 


s . ) Re 

l 



(k = 1, 2, . . . , M) 


(71) 


Further simplication is possible if s is a memory less mapping; 
thus, with x^ = x(i/2B), Eq. (69) becomes 


V x £nA 


2a 


1 + aN 


r * 9s i" 

(k . , s . ) Re s . 7T 

i' i l 9x. 

L 1 J 


( i = 1, 2, 


2Bt) 


(72) 


In general, it appears that a search procedure 


** — j 


find x.,,._. or x.. T . 
MAP ML 


If we now assume that the a priori density P x (^) [x(t)] is 
gaussian, then Eq. (67) reduces to 
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Notice that the optimum MAP estimate is performed on the intensity 
of the process. That is k^ is first compared to a | s.jj , the 
expected count of the received process. Since | | is the 
envelope of the process, it can be written as e^ n l s il, and 
gT £n j s^ | serves the same function as yg- s(t, x) in the MAP 
equations for gaussian systems (ref. 19, p. 432). 

If the intensity of the process is constant, such as FM, 

3 

no phase information can be obtained since 3 — In | s . | = 0 , making 

o X 1 

x = 0. To obtain phase information from an FM or a PM signal, 
optimal heterodyne detection must be employed. That is, the 
signal from a local oscillator is aligned and mixed with the 
received signal over the surface of the detector. Then, 


s . 
1 


E. 


, ju), t + <p (t) 


+ E, 


3“ 2 t 


where and E Lq are the two electric field strengths, and 





o 


* D ( 

E t e 
-Li 

o 

- w 2 ) 


W 1 " w 2 


t + <J>] 


'} 


) t + <f>J 


(76) 
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Assuming narrow-band signals, we have 


s.| 2 - | E | 2 + |E I 2 + 2E E cos 

o o 


By making the local oscillator signal large | E L j >> ]e^ 


| 2 £l + g— ^ cos - w 2 )t + <f>(t)jj 


1 2 

Zn | s . | = 2 Zn | s^ | - £n |e | + - cos 

1 Z 1 o 

o 



We see, then, that the optimum detector is similar, in this case, 
to a phase-lock loop. It can also be shown by using Eq . (26) and 
references 18 and 26 that if we consider the random variable k^ 
as a shot-noise process and pass it through a narrow-band filter 
with bandwidth 2W tuned to (co-^ - w 2 ) / the density will approach 
a gaussian density with mean 2aE^E^ cos [ (u)^ - w 2 )t + 4> ( t ) ] and 
variance a|E LQ |2W. 

Using the notion of an information kernel , we can lower- 
bound the mean-square error of any estimator resulting from Eq. 

(66) , without specifying the form of the estimator (ref. 19, p. 80) . 
Defining 


- th P^ Zn__A 3 £n A ~| 
J D - |^3x(u) 3x(v) J 


(80) 
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( 81 ) 


= E 


C 


’3 In p 


x (t ) 


[X (x) ] 3 £n p 


3x (u) 


x (t) 


[X(t)] 


3x (v) 


] 


J T = J D + Jp = information kernel 
we can show that 


( 82 ) 



[x(x) 


x (t) ] 2 dx ^ Tr J T 1 


( 83 ) 


The expectation in Eqs. (80) and (81) is over the random function, 
x. Combining Eqs. (69) and (80), we have 


J D = E« 


4a" 


(1 + aN Q ) i=l j=l 


2Bt 2Bt 

2 E E = <V S.) 


r * 9s i 1 

f * 9s i 1 

i 3x(u) 

L j 3x (v) 


( 84 ) 


To evaluate this further, the transformation s must be specified. 


CONCLUSIONS 


We have presented general results for the photoelectron 
counting distributions arising from the quantum detection of a 
narrow-band gaussian process. The semiclassical approach has led 
to probability distributions which could be written explicitly 
in terms of the covariance function of the gaussian radiation. 
Four different, but equivalent, methods for evaluating the 
counting distribution have been presented, the last of which is 
valid only for white noise. These could be descriptively titled, 
"the eigenvalue approach," "the iterated kernel approach," "the 
compound Poisson approach," and "the time-sampling approach." 

The third method is the one in which Eq. (3) is evaluated 
directly. Of the four methods, the third and last one lend the 
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most physical insight, although the third method is seldom tracta- 
ble, computationally. The eigenvalue approach will probably find 
the most use in problems of practical importance. 

Several examples have been given, with special emphasis on 
the important case of a deterministic signal in white gaussian 
noise. Some limiting forms have been derived, and rigorous 
justification has been given for the often-made assumption of a 
Poisson signal in Poisson noise. Finally, the likelihood func- 
tional has been defined, and formal answers have been found for 
the continuous MAP and ML estimation problems, including a simple 
approximation for low noise, and an expression for the minimum 
mean-square error of any estimator. 

It should be emphasized that, from a communications view- 
point, the detector senses variations in the signal modulus only; 
all carrier phase information is lost. This does not mean, 
however, that our results are applicable only to direct (or 
"energy") detectors. Heterodyne detection can be accomplished 
simply by illuminating the detector with a local oscillator field; 
the squared modulus then contains the dif ference frequency term 
and the phase information from both the local oscillator and the 
received field. 
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APPENDIX 


Here we mention a few of the more important statistical 
quantities associated with the counting distribution, p N ^ (k) . 
First we note that the factorial moments ^ [ jj of the counting 

distribution are equal to the moments {y^} of m t . This is clear 
from 


y r , = X) k(k - 1) . . . (k - n + l)p (k) 
1 J k=n w t 


n 9 


n 


= (-1)“ M [Jin (1-5)1 

95 W t 


K = 0 


,n 


(_l) n 1 M (~5) 

n m 

9? t 


5=0 


= u n • (Ai) 

Also, the cumulants of the counting distribution are connected 
with the cumulants of m by 

n k . 

< = L A (n, i) -A- (A2) 

n i=l 1 * 


where 


A (n , i) 



(A3) 


*The tilde (~) will be used to distinguish moments, etc., of the 
counting distribution, from corresponding quantities associated 
with m t * 
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Equation (A2) is a consequence of differentiating the logarithm 

of 


(jv) = M (e jv - 1) 
t t 


= exp E (e^ V - l) k 

k=l 


(A4) 


Similarly, it can be shown that 


y 


n 


n y . 

E A (n , i) 
i=l 1 * 


(A5) 


Further, the moments {y^} of are related to the cumulants 

{k^} of by Eq. (35) . 

With the help of these formulas and known identities (ref. 
24) , an exhaustive set of relations can be found for the moments, 
factorial moments, cumulants, and factorial cumulants of and 
m^. We shall not, however, reproduce these relations here. 
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